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摘要 : 青冈 亚 属 植物 的 系统 发 育 地 位 长 期 存在 着 争议 ， 部 分 种 的 种 间 关 系 不 明确 。 为 揭示 宁 
内 青冈 (Quercus ningangensis) 、 曼 青冈 (QO. oxyodon) 、 毛 曼 青冈 〈O. gambleana? ~ 4y 
ERK CQ. neglecta) 的 叶绿体 基因 组 特征 及 系统 发 育 关 系 ， 本 研究 选择 以 上 4 种 青冈 的 成 
熟 叶片 进行 二 代 测 序 , 对 其 叶绿体 基因 组 结构 和 特征 进行 分 析 , 并 结合 相关 类 群 进行 系统 发 
AWA. ARRA: d) 宁 冈 青 网 、 曼 青冈、 毛 曼 青冈 、 竹 叶 青 冈 的 叶绿体 基因 组 序列 长 
度 分 别 为 : 160 906 bp. 160 883 bp、160 832 bp. 160 784 bp， 它 们 均 编码 133 个 基因 ， 包 括 
88 个 蛋白 质 编 码 基因 、37 个 tRNA 基因 、8 个 rRNA 基因 。 〈2) 四 种 青冈 更 偏好 以 AT 结 
尾 的 密码 子 ， 质 体 基 因 组 变异 区 域 主要 存在 于 非 编码 区 域 。 (3) 通过 IR 边界 分 析 得 出 ， 四 
PE pq NV S SLE TE ycfl 假 基 因 ， 且 在 IRb/SSC 区 域 发 生 扩张 。 (4) 系统 发 育 分 析 表 明 ， 
ETAPER, 水 青冈 属 (Fagus) 和 三 棱 栎 属 (Trigonobalanus) 较 早 分 化 出 来 , 栎 亚 属 Csubg. 
Quercus) 未 形成 一 个 单 系 群 ， 叶 绿 体 基因 组 建树 结果 和 核 标 记 、 质 体 标记 一 致 ， 栎 亚 属 中 
Ilex 组 和 Cerris 组 间 出 现 穿插 现象 。 (5) 《中 国 植物 志 》 及 部 分 地 方 植物 志 将 毛 曼 青 闷 作 
为 一 个 单独 的 种 处 理 ， 而 周 浙 昆 、 邓 敏 等 人 根据 叶 的 特征 将 毛 曼 青冈 作为 曼 青 冈 的 一 个 亚 种 
处 理 ， 毛 曼 青 冈 的 系统 发 育 地 位 存在 着 争议 ,本 实验 根据 叶绿体 基因 组 信息 并 结合 前 人 的 形 
态 分 析 结 果 , 支持 毛 曼 青冈 作为 一 个 独立 的 种 存在 的 观点 。 该 研究 为 探讨 青冈 亚 属 系统 发 育 
地 人 位、 青冈 亚 属 组 的 划分 、 存 疑 种 种 间 关 系 的 解决 提供 了 基础 资料 。 
关键 词 : 老 斗 科 ， 青 冈 亚 属 ， 基 因 组 比较 ， 叶 绿 体 基因 组 ， 系 统 发 育 研 究 


Comparison of chloroplast genomes and phylogenic analysis 


of four Species of Quercus subg. Cyclobalanopsis 
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Abstract: The phylogenetic status of Quercus subg. Cyclobalanopsis has long been controversial 
and the interspecific relationships of some species are unclear. At present, chloroplast genomes 
have been used to solve the phylogenetic problems of subg. Quercus and Castanea in the 
Fagaceae, but there are few researches on solving the phylogenetic problems of Quercus subg. 
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Cyclobalanopsis with chloroplast genome information. To reveal the chloroplasts genome 
characteristics and phylogenetic relationships of Q. ningangensis, Q. oxyodon, Q. gambleana and 
Quercus neglecta. In this study, mature leaves of the above four species of Quercus subg. 
Cyclobalanopsis were selected for Next-generation sequencing, and the chloroplast genome 
structure and characteristics were analyzed, Phylogeny was studied in combination with related 
taxa. The results were as follows: (1) The chloroplast genome sequences of Q. ningangensis, Q. 
oxyodon, Q. gambleana and Q. neglecta were 160 906 bp, 160 883 bp, 160 832 bp and 160 784 
bp, respectively. They all encoded 133 genes, including 88 protein-coding genes, 37 tRNA genes 
and 8 rRNA genes. (2) The codons ending in A/T were preferred by the four species of Quercus 
subg. Cyclobalanopsis, and the variation regions of plastid genome mainly existed in non-coding 
regions. (3) According to IR boundary analysis, yc/f7 pseudogene was found in four species of 
Quercus subg. Cyclobalanopsis and expanded in IRb/SSC regions. (4) Phylogenetic analysis 
showed that Fagus and Trigonobalanus were differentiated early in Fagaceae, while subg Quercus 
did not form a monophyly. The results of phylogenetic tree based on chloroplast genome were 
consistent with the nuclear marker plasmid markers. Interspersed between Ilex group and Cerris 
group in subg Quercus. (5) The Flora of China and some local flora treated Q. gambleana as a 
separate species, while Zhou Zhekun and Deng Min et al. treated Q. gambleana as a subspecies 
according to the characteristics of leaves. The phylogenetic status of Q. gambleana is still 
controversial. Based on the chloroplast genome information and previous morphological analysis 
results, the present study supported the idea that Q. gambleana existed as an independent species. 
This study provides basic data for discussing the phylogenetic status of the Quercus subg. 
Cyclobalanopsis, the division of its groups, and the resolution of doubtful relationships among 
them. 
Keywords: Fagaceae, Quercus subg. Cyclobalanopsis, genome comparison, chloroplast genome, 
phylogenic analysis 

叶绿体 广泛 存 于 高 等 植物 及 大 部 分 藻类 中 , 是 绿色 植物 进行 光合 作用 的 场所 。 叶绿体 自 
身 拥有 独立 的 遗传 物质 ， 其 基因 组 为 四 分 体 结构 ， 由 一 个 大 单 拷贝 区 域 Cong single copy, 
LSC) .一 个 小 单 描 贝 区 域 (small simple copy, SSC) 及 两 个 反 向 重复 区 域 (inverted repeat, IRA 
and IRB) 构成 ， 基 因 组 序列 和 结构 都 相对 保守 (Zhao et al., 2018) 。 与 核 基因 组 比较 结果 
显示 ,叶绿体 基因 组 具有 结构 简单 、 分 子 量 小 、 基 因 的 替换 率 较 低 、 单 拷贝 等 特点 (Korpelainen, 
2004; Kwak et al., 2019) 。 近 年 ， 随 着 二 代 测 序 时 间 和 成 本 的 降低 ， 越 来 越 多 的 叶绿体 基因 
昌 的 测序 分 析 结 果 已 经 用 于 模糊 种 的 鉴定 和 系统 发 育 研 究 。 目 前 ， 已 有 研究 利用 叶绿体 基因 
日 解决 壳 斗 科 栎 亚 属 Csubg. Quercus? ~ RJE (Castanea) WIZ&Z AB WE CEA, 2018; 
高 潇潇 , 2020) ， 但 利用 叶绿体 基因 组 信息 解决 青冈 亚 属 系统 发 育 问题 的 还 未 多 见 。 

壳 斗 科 (Fagaceae) He E KEE (Quercus subg. Cyclobalanopsis) 共有 122 种 (Frodin 
& Govaerts, 1998) 或 150 种 (Huang et al., 1999) ， 主 要 分 布 于 亚洲 热带 和 亚热带 地 区 。 我 
国 是 青冈 亚 属 的 重要 分 布 区 和 多 样 性 中 心 之 一 〈 罗 艳 和 周 浙 昆 , 20010 ， 共 有 77 种 以 及 3 
变种 ， 秦 岭 、 淮 河 以 南 为 青冈 亚 属 植物 的 主要 分 布 区 ， 为 常 绿 阔 叶 林 的 建 群 种 之 一 〈 陈 焕 饼 
和 黄 成 就 , 1998; Huang et al., 1999) 。 青 办 类 植物 木材 具有 质地 坚硬 、 耐 腐蚀 、 耐 磨损 等 优 
点 ， 是 制作 船 奖 和 其 他 工具 的 良好 木材 ， 其 果实 内 淀粉 含量 丰富 ， 含 量 约 为 5096-6096 Cin 
木 煌 , 1995) ， 可 用 于 酿酒 、 制 作 饲 料 、 糊 料 和 工业 淀粉 等 ， 树 皮 中 含有 单 守 ， 可 制作 楼 胶 。 
青冈 类 植物 的 系统 发 育 地 位 长 期 存在 着 争议 , 部 分 学 者 根据 壳 斗 的 排列 方式 将 青冈 类 植物 作 
为 一 个 单独 的 属 处 理 〈 徐 永 椿 和 任 宪 威 , 1976; 郑 万 钧 , 1985; 陈 焕 镭 和 黄 成 就 , 1998; Huang 
et al., 1999) ; 但 多 数学 者 根据 花粉 粒 解剖 学 等 形态 特征 及 ITS 序列， 支持 青 冈 类 植物 作为 
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栎 属 下 的 一 个 亚 属 处 理 ( 中 国 科 学 院 植物 研究 所 , 1972; 刘 兰 芳和 房 志 坚 , 1986; 端木 煌 , 1992; 
罗 艳 和 周 浙 昆 , 2001b) 。 由 于 青冈 类 植物 生殖 结构 变化 较 大 ， 叶 的 特征 较为 均一 ， 不 同类 群 
间 存 在 基因 渐 渗 ， 和 馆藏 的 标本 量 有 限 ， 果 实 等 重要 分 类 特征 难以 观察 ( 邓 敏 ,2007) ， 故 其 
盟 下 没有 明确 的 分 组 ， 部 分 种 的 种 间 关 系 不 明确 ,种 的 划分 上 还 存在 着 许多 争议 。 目 前 关于 
青冈 亚 属 系统 发 育 的 研究 主要 还 停留 在 形态 学 方面 。 因此 , 需要 更 多 的 分 子 系统 学 方面 的 证 
据 来 协助 解决 青冈 亚 属 的 系统 发 育 地 位 和 种 间 关 系 。 

本 研究 对 4 种 青冈 亚 属 植物 的 成 熟 叶片 进行 二 代 测 序 , 依托 4 种 青冈 亚 属 植物 的 叶绿体 
基因 组 信息 ， 拟 探究 以 下 问题 : 〈1) 4 种 青冈 亚 属 植物 的 叶绿体 基因 组 基本 特征 及 其 差异 。 
(2) 4 种 青冈 亚 属 植物 的 密码 子 使 用 偏好 性 。 (3) 毛 曼 青冈 和 曼 青冈 亲缘 关系 如 何 ? 以 期 
为 青冈 亚 属 系统 发 育 问题 解决 、 种 间 关 系 的 划分 、 种 的 鉴定 及 优良 青冈 木材 的 选 育 从 叶绿体 
基因 组 层面 提供 了 新 的 资料 。 


1 材料 与 方法 
1.1 实验 材料 


EEI RETAK (Quercus ningangensis) 、 曼 青冈 CQ. oxyodon) ~- ESF] O. 
gambleana) . TW ERA] CQ. neglecta) 新 鲜 的 、 无 病虫害 的 成 熟 叶 片 〈 表 1) ， 用 湿润 的 纱 
布 擦拭 后 置 于 装 有 变色 硅胶 的 密封 袋 中 保存 ， 干 燥 后 置 于 -80C 的 冰箱 内 备用 。 和 凭证 标本 存 
放 于 西 华师 范 大 学 标本 馆 CHerbarium, College of Life Sciences, China West Normal 
University) 。 
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表 14 种 青冈 亚 属 植物 样品 采集 信息 
Table 1 Information on the collection of samples from four species of 


Quercus subg. Cyclobalanopsis 


植物 名 采集 地 点 经 纬度 海拔 (m) 凭证 标本 号 
Plant name Collection place Latitude/Longitude Altitude (m) Voucher specimen 
TAEK T 
WARTH N26°22'24.74" 
Quercus 1233 HT211016 
Suining, Hunan E110?07'54.40" 


ningangensis 


z 广西 兴安 县 猫 儿 山 
曼 青 区 N25?5421.67" 


Maoer Mountain, 149 8 HT210925 
Q. oxyodon . . E110?27'58.79" 
Xingan, Guangxi 
毛 曼 青冈 重庆 市 南川 区 N29*00'58.46" 
T HT211014 
. gambleana anchuan, Chongqin ; 
gambl Nanch Chongqing E107°11'18.90” 
E I: 全- vh E v4 
SUN C EE L 
竹叶青 冈 i N27953'48.29" 
Fanjin Mountain, 1563 HT210917 
Q. neglecta E108?43'19.03" 


Tongren, Guizhou 

1.2 基因 组 DNA 的 提取 与 测序 

用 改良 的 CTAB 法 从 4 种 青冈 亚 属 植物 叶片 中 提取 总 DNA。 利 用 Humina Hiseq 对 提取 
的 DNA 进行 双 端 测序 。 叶 绿 体 全 基因 组 DNA 的 提取 与 测序 均 由 南京 奥 维 森 基因 科技 有 限 
公司 完成 。 
1.3 组 装 与 注释 

利用 FASTQ 软件 进行 数据 质 控 , 将 质 控 得 到 的 clean data 利用 GetOrganelle 软件 (Jin et 
al., 2020) 进行 组 装 ， 将 原始 数据 map 到 组 装 好 的 序列 上 查看 履 盖 率 并 结合 Gapcloser 软件 
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(Zuo etal., 2017) ， 补 齐 组 装 出 现 的 空缺 (gaps) o Æ NCBI (https:/www.ncbi.nlm.nih.gov/) 
E PAS PX 《MN199023) 的 叶绿体 基因 组 作为 参考 基因 组 ， 利 用 CPGAVAS2 (Shi et al., 
2019) (http://www.herbalgenomics.org/cpgavas/ ) 进行 注释 ， 注 释 好 的 数据 利用 Geneious 软 
件 C Kearse et al, 2012 ) 进行 手动 调整 和 修改 。 利 用 OGDRAW 软件 
(nttps://chlorobox.mpimp-golm.mpg.de/OGDraw.html) 绘制 4 种 青冈 亚 属 植物 的 叶绿体 基因 
晶 物 理 结构 图 谱 。 最 后 将 注释 好 的 数据 提交 到 NCBI TAEK SFK EZEK r 
和 四 的 序列 登录 号 分 别 为 : ON303301、ON258628、ON258629、ON258631。 

1.4 叶绿体 基因 组 比较 分 析 
利用 Geneious (Kearse et al., 2012) 及 Editseq 软件 统计 4 种 青 内 亚 属 植物 的 叶绿体 基因 
组 序列 SSC、LSC、IRS 区 域 的 长 度 、GC 含量 等 基本 信息 。 从 NCBI 上 下 载 小 叶 青冈 
(MN199025) 、 神 时 青冈 (MN199023 ) . Pu (NC 036941) 、 青 冈 (NC _036930) 
的 叶绿体 基因 组 与 我 们 测 得 四 种 青冈 亚 属 植物 的 叶绿体 基因 组 进行 比较 分 析 。 
选择 以 ATG 为 起 始 密码 子 、 长 度 大 于 300 pb 的 非 重 复 序列 ， 所 有 物种 筛选 后 均 剩 余 52 
条 CDS 序列 。 利 用 Codon W 软件 (Sharp & Li, 1987) 计算 4 种 青冈 亚 属 植物 及 其 近 缘 种 的 
同 义 密码 子 相 对 使 用 度 Crelative synonymous codon usage, RSCUO 、 密 码 子 适应 指数 (codon 
adaptation index, CAD 、 有 效 密码 子 数 (effective number of codons, ENC) 、 密 码 子 偏好 性 
指数 (codon bias index, CBI) 、 最 优 密码 子 使 用 频率 (frequency of optimal codons, FOP) 以 
及 每 条 CDS 序列 的 ENC 值 、RSCU 值 及 密码 子 第 三 位 各 碱 基 出 现 的 概率 〈 记 为 A3、T3、 


uL Me 


C3. G3) ， 运 用 在 线 软件 EMBOSS ( https://www.bioinformatics.nl/emboss-explorer/ ) 计算 各 


CDS 序列 GC1、GC2、GC3 的 含量 。 中 性 绘图 分 析 以 GC1 和 GC2 的 平均 值 CGC) 为 纵 
坐标 ， 以 GC3 值 为 横 坐 标 ， 绘 制 散 点 图 ， 并 插入 yox 的 函数 图 像 。ENC-plot 分 析 绘 制 的 散 
点 图 将 GC3 值 、ENC 值 分 别 作 为 横 、 纵 坐标 ， 添 加 理论 ENC 值 曲线 ， 标 准 曲 线 公 式 为 
ENC-24GC3429/[GC3^4(1-GC3y]. PR2-plot 分 析 以 G3/(G3+C3) 为 横 坐 标 ，A3/(A3+T3) 为 纵 
坐标 绘制 散 点 图 , 分 析 A、T、C、G 的 使 用 频率 和 偏 倚 性 。 中 性 绘图 、ENC-plot 22 Ed . PR2-plot 
绘图 均 使 用 origin 软件 。 

利用 mVISTA 软件 (https://genome.lbl.gov/vista/mvista/instructions.shtml) (Mayor et al., 
20000 对 4 种 青 闷 亚 属 植物 及 其 近 缘 种 的 叶绿体 基因 组 进行 可 视 化 对 比分 析 。 使 用 IRscope 

Chttps://irscope.shinyapps.io/irapp/) (Amiryousefi et al., 2018) AFREK SFK E 

曼 青冈 、 竹 叶 青 冈 及 其 近 缘 种 的 叶绿体 基因 组 IR 边界 。 
1.5 系统 发 育 分 析 

选取 并 下 载 在 NCBI 中 公开 释放 的 8 种 壳 斗 科 植 物 的 叶绿体 基因 组 序列 , 包括 : 倒 卵 叶 
青冈 (QO. obovatifolia, MG356785.1) 、 川 溃 高 山 标 《OQ. aquifolioides, KX911971.1) 、 刺 叶 高 
Ii#k CO. spinosa, MG678038.1) 、 巴 东 栎 (QO. engleriana, MZ196209.1) ~ {EEk CQ. variabilis, 
NC031356.1) ~ Z}REk (Trigonobalanus doichangensis, NC0239591) 、 人 台湾 水 青冈 、 米 心 水 
青冈 选用 台湾 水 青冈 Fagus hayatae) 和 米 心 水 青冈 CE engleriana) 作为 外 类 群 ， 将 上 述 8 
种 壳 斗 科 植 物 的 叶绿体 基因 组 与 本 文 组 装 的 4 种 青冈 亚 属 植 物 的 叶绿体 基因 组 一 起 构建 系 
统 发 育 树 。 将 所 有 叶绿体 基因 组 序列 运用 MAFFT (Katoh & Standley, 2013) 进行 对 比 ， 将 
对 比 好 的 叶绿体 基因 组 序列 利用 MEGA 7 软件 (Kumar et al., 2016) 对 序列 开头 和 结尾 未 对 
齐 部 分 进行 适当 的 修剪 ， 运 用 MEGA 7 软件 中 Models 功能 寻找 以 上 序列 运用 ML 法 构建 系 
统 发 育 树 的 最 佳 模 型 (GTR+GI) ,以 最 大 似 然 法 Gnaximum likelihood, ML) ,设置 bootstrap 
为 100 0 构建 系统 发 育 树 。 另 将 修剪 好 的 序列 再 运用 MEGA 7 软件 以 邻接 法 Cneighbor-joining 
method, NJ) ， 选 用 Maximum Composite Likelihood 模型 ， 设 置 bootstrap 为 100 00 构建 相关 
物种 的 系统 发 育 树 。 
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2 结果 与 分 析 


2.1 叶绿体 基因 组 的 基本 特性 


TAEK ZF K, 62 


青冈 、 竹 叶 青 冈 的 叶绿体 基因 组 均 有 1 


ChinaXiv 合 作 期 刊 


一 个 LSC 


区 域 、 两 个 


IRs 区 域 、 一 个 SSC 区 域 构 成 , 为 典型 的 四 分 体 结 构 (图 10 。 其 基因 组 长 度 为 160 681~160 


906 bp， 其 中 曼 青 内 和 毛 曼 青冈 叶绿体 全 基因 


组 的 长 度 分 别 为 160 883 bp 和 160 882 bp. fX. 


相差 1 个 碱 基 。LSC 区 的 长 度 范围 为 90 245-90 360 bp, SSC 区 的 长 度 范围 为 18 891-18 929 
bp，IRs 的 长 度 范围 为 25 816-25 840 bp。 本 实验 测 得 的 四 种 青冈 的 总 GC 含量 (36.9%) 、 


IRs 区 域 的 GC 含量 (42.8%) 及 CDs 区 域 的 GC 含量 
区 域 的 GC 含量 也 相差 不 大 R 2) 。 


(37.9%) 均一 致 ，LSC 区 域 及 SSC 


宁 冈 青冈 、 曼 青冈 、 毛 曼 青冈 、 竹 叶 青 冈 的 叶绿体 基因 组 注释 结果 显示 ，4 REDDE JS 


植物 均 有 133 个 基因 ， 其 中 包含 37 个 tRNA 基因 Ctransfer RNA genes) ，8 个 rRNA 基因 
(ribosomal RNA genes) ， 以 及 88 个 和 蛋白质 编码 基因 CProtein-coding Genes? ， 
有 15 个 基因 有 1 个 内 售 子 GpoCI. ndh4、 ndhB、 rpl2. rpl16、 atpF、 rps16、 trnL-UAA, 


( 表 2) 。 


trnK-UUU, trnl-GAU, trnA-UGC, trnG-GCC, trnV-UAC, petB, petD) ， 包 含 2 个 内 含 子 
yeso ( 表 3)。 


的 基因 有 3 个 (rps12、clpP、 


PE 


圆圈 内 外 的 灰色 箭 
灰色 表示 AT 含量 。 


Quercus ningangensis 160 906 bp 
曼 青 冈 
Q. oxyodon 160 883 bp 
毛 曼 青冈 
Q. gambleana 160 832 bp 
竹叶青 冈 
Q. neglecta 160 784 bp 


E DA] 


色 表 示 GC 含量 1€ 


Gray arrows inside and outside the circle indicate the direction of gene transcription; Different colors represent 


genes with different functions; Dark gray in the inner circle represents GC content, the light gray represents AT 


content. 


图 


1 4 种 青冈 亚 属 植物 的 叶绿体 基因 


组 图 


hi 


it 


H 


Fig. 1 Gene map of chloroplast genome in four species of Quercus subg. Cyclobalanopsis 


de 24 种 青冈 亚 属 植物 叶绿体 基因 组 特征 比较 
Table 2 Comparison of chloroplast genomes among four species of 


Quercus subg. Cyclobalanopsis 


表 3 4 fires Pu E Je fW Ee s ks DSL 2H S DSL 2L 


物种 FAEN PETS 曼 青冈。 竹叶青 网 
Species Quercus ningangensis Q. oxyodon Q. gambleana Q. neglecta 
基因 组 大 小 
160 906 160 883 160 832 160 784 
Genome size (bp) 
大 单 拷贝 区 长 度 
90 297 90 360 90 259 90 245 
Length of LSC (bp) 
小 单 拷贝 区 长 度 
18 929 18 891 18 905 18 895 
Length of SSC (bp) 
反 向 重复 区 长 度 
25 840 25 816 25 834 25 822 
Lengh of IRs (bp) 
蛋白 编码 基因 数 
: 88 88 88 88 
Number of protein-coding genes 
转运 RNA 基因 数 
37 37 37 37 
Number of tRNA genes 
核糖 体 RNA 基因 数 à i " " 
Number of rRNA genes 
总 基因 数 
133 133 133 133 
Total number of genes 
总 的 GC 含量 
36.9 36.9 36.9 36.9 
GC content of overall (96) 
AGERE WIX GC 含量 
34.8 34.8 34.7 34.8 
GC content of LSC (96) 
小 单 拷贝 区 GC 含量 
31.1 31.1 31.1 31.0 
GC content of SSC (96) 
反 癌 重复 区 GC 含量 
42.8 42.8 42.8 42.8 
GC content of IRs (%) 
基因 编码 区 GC 含量 
37.9 37.9 37.9 37.9 


GC content of CDs (%) 
Table 3 Composition of complete chloroplast genome of four species of 
Quercus subg. Cyclobalanopsis 


基因 功能 分 类 EE [A] 4] 2H. 基因 名 称 数量 
Category of gene Group of gene Name of gene Number 
Te rrn4. 5S( x2); rrn58(*2); rrnl6S(x 
自我 复制 核糖 体 RNA EAR : 
Self-replication rRNA genes 


2); rrn23S(x2) 


光合 作用 
Photosynthesis 


转运 RNA 基因 
tRNA genes 


RNA JE 


RNA polymerase 


核糖 体 小 亚 基 


Small subunit of ribosomal 


核糖 体 大 亚 基 


Large Subunit of ribosomal 


NADH- 脱 氧 酶 
NADH dehydrogenase 


光 系 统 1 


Photosystem I 


光 系 统 2 
Photosystem II 


细胞 色素 bf 复合 物 


trnA-UGC(X2)*; trnC-GCA; 


trnD-GUC; trnE-UUC; trnF-GAA; 
trnfM-CAU; trnG-GCC*; 
trnG-UCC; trnH-GUG; 


trnl-CAU(X2); trnl-GAU(X2)*; 


trnK-UUU*; trnL-CAA( X2); 


trnL-UAA*; trnL-UAG; 


trnM-CAU; trnN-GUU( x2); 
trnP-UGG; trnQ-UUG; 


trnR-ACG(X2); trnR-UCU; 


trnS- GCU;trnS-GGA; 


trnS-UGA;trnT-GGU; trnT-UGU; 


trnV-GAC(*2); trnV-UAC*; 


trnW-CCA; trnY-GUA 


rpoA; rpoB; rpoC1*; rpoC2 


TDS2; rps3; rps4; rps7(X2); rps8; 


rpsll; rps12(x2)**; 


rps14; rps15; rpsl6*; rps18; rps19 


rpl2(x2)*; rpl14; rpl16*; rpl20; 


rpl22; rpl23(x2); 


rpl32; rpl33; rpl36 


ndhA*; ndhB(x2)*; ndhC; ndhD; 


ndhE; ndhF; 
ndhG; ndhH; ndhlI; ndhJ; 
ndhK 


psaA; psaB; psaC; psal; psaJ 


psbA; psbB; psbC; psbD; psbE; 


psbF; psbH; psbl; 
psbJ; psbK; psbL; psbM; 
psbN; psbT; psbZ 


petA; petB*; petD*; petG; petL; 


37 


14 


11 


12 


15 


Cytochrome b/f complex petN 


ATP A FX atpA; atpB; atpE; atpF*; atpH; P 
ATP synthase atpl 
二 磷酸 核 酮 糖 羧 化 酶 大 亚 基 基因 T i 
rbc. 
Large subunit of rubisco 
蛋白 酶 基因 
clpP** 1 
Protease 
成 熟 酶 基 攻 
matK 1 
Maturase 
乙酰 CoA 羧 化 酶 的 亚 基 
Subunits of accD 1 
其 他 基因 
Acetyl-CoA-carboxylase 
Other genes 2 
包 膜 蛋白 基因 
cemA 1 
Envelope membrane protein 
C 型 细胞 色素 合成 
. ccsA 1 
C-type cytochromesynthesis 
翻译 起 始 因子 
mE infA 1 
Translation initiation factor 
未 知 功能 基因 D ycf1(X2); ycf2(x2); ycf3**; ycf; 
开放 阅读 框 
Unknown function 8 
Open reading frames 
gene ycf15(x2) 
总 数 Total 133 


注 : (X2) 两 次 拷贝 ，* 一 个 内 含 子 ; ** 两 个 内 含 子 。 
Note: (x2) two copies; * One intron; ** Two introns. 


2.2 密码 子 偏好 性 分 析 


2.2.1 密码 子 组 成 分 析 


4 种 青冈 亚 属 植物 的 密码 子 数 在 20 996~20 977 ZE; 有 效 密码 子 数 (effective number of 
codons, ENC ) 在 49.91~49.81 之 间 , 密码 子 适 应 指数 (codon adaptation index, CAT) 均 为 0.167; 
宁 冈 青冈、 毛 曼 青冈 和 竹叶青 网 的 密码 子 偏好 性 指数 (codon bias index, CBI) 73-0.099, 6 
青 内 的 密码 子 偏好 性 指数 为 -0.097; 4 种 青冈 亚 属 植物 的 最 优 密码 子 使 用 频率 (frequency of 
optimal codons, FOP) 均 为 0.355，GC 含量 为 37.93%~37.95% 之 间 。4 种 青冈 亚 属 植物 的 密 
码 子 偏好 性 各 指数 相差 不 大 ， 说 明 它 们 密码 子 使 用 情况 相似 〈 表 4) 。4 种 青冈 亚 属 植物 的 
GC3 含量 在 29.85%~29.88% 之 间 ， 表 明 青 闷 亚 属 植物 更 偏好 以 AT 结尾 的 密码 子 。 

如 图 二 , ZAR (Leu)、 丝 氨 酸 Ser) 精 氨 酸 Arg) 有 6 个 同 义 密码 子 , CAR (Trp)、 
HAE? (Met) 仅 由 一 种 密码 子 编码 ， 其 余 氨 基 酸 均 由 2 个 及 以 上 的 同 义 密码 子 编码 。 这 
8 种 青冈 亚 属 植物 的 密码 子 偏 好 性 相似 ， 其 中 ACU、UCU、UUA、GCU、UAU、GAU、 

AGA 的 使 用 频率 较 高 , CUC, CUG, GCG, UAC, CAC, CAG, AAC, AAG, GAC, GAG, 
CGC, CGG, AGC, GGC 的 使 用 频率 较 低 。 在 米 楼 (Castanopsis carlesii) , KRA (Lithocarpus 
longinux) 中 UAA 的 使 用 频率 较 所 列 的 其 他 植物 低 。 
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严 属 植物 的 密码 子 偏好 性 指数 


Table 4 Codon preference index of four species of Quercus subg. Cyclobalanopsis 


指数 


Index 


长 度 
Length (bp) 
密码 子 数 
Codon number 
有 效 密码 子 数 
Effective number of codons 
密码 子 适 应 指数 
Codon adaptation index 
密码 子 偏好 指数 
Codon bias index 
最 优 密码 子 使 用 频率 
Frequency of optimal codons 
GC 含量 
GC content (%) 
GCI 含量 
GCI content (%) 
GC2 含量 
GC2 content (96) 
GC3 含量 


GC3 content (96) 


TAEK 


Quercus ningangensis 


63 021 


20 990 


49.89 


0.167 


-0.099 


0.355 


37.93 


46.12 


37.84 


29.85 


Q. 


注 : 白色 白色 和 红色 分 别 表示 较 低 和 


较 高 的 RSCU E. PSEK SAK Yr 
组 装 和 注释 ， 其 余 物 种 均 为 NCBI 上 下 载 。 


Note: White and red indicate lower and higher RSCU values, respectively. Quercus gambleana, Q. oxyodon, 


曼 青冈 毛 曼 青冈 竹叶青 冈 
oxyodon | Q.gambleana Q. neglecta 
62 982 62 991 63 039 
20 977 20 980 20 996 
49.90 49.91 49.9] 
0.167 0.167 0.167 
-0.097 -0.099 -0.099 
0.355 0.355 0.355 
37.95 37.94 37.95 
46.12 46.11 46.12 
37.87 37.85 37.85 
29.87 29.87 29.88 


台湾 水 青冈 Fagus hayatae MW846258.1 
光 叶 水 青冈 F. lucida NC043858.1 
米 心 水 青冈 F. engleriana NC036929.1 

三 楼 栎 Trigonobalanus doichangensis NC023959.1 
FÆ Castanea seguinii MK248944.1 

西 畴 青 网 Quercus sichourensis NC036941.1 
青冈 O. glauca NC036930.1 

毛 曼 青 网 O. gambleana ON258629 
SHP] Q. oxyodon ON258628 

答 叶 青冈 O. neglecta ON258631 

褐 叶 青冈 Q. stewardiana MN199023.1 

宁 冈 青冈 Q. ningangensis ON303301 

小 叶 青冈 Q. myrsinifolia MN199025.1 

栓 皮 栎 Q. variabilis NC031356.1 

巴 东 榴 Q. engleriana MZ196209.1 

MIRI Q. dentata NC039725.1 

蒙古 栎 O. mongolica NC043858.1 

米 村 Castanopsis carlesii MK745999.1 

长 果 柯 Lithocarpus longinux OK181903.1 


内 、 宁 冈 青 冈 为 该 研究 


Q. neglecta and Q. ningangensis were assembled and annotated for the study, and the remaining species were 


downloaded from NCBI. 


图 2 4 种 青冈 亚 属 植物 及 近 缘 种 的 同 义 密码 子 的 相对 使 用 度 
Fig.2 RSCU values of all protein-coding (CDS) genes for chloroplast genomes from four 


species of Quercus subg. Cyclobalanopsis and its sibling species 
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2.2.2 中 性 绘图 分 析 
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从 4 种 青冈 亚 属 植物 叶绿体 基因 组 筛选 出 以 ATG 为 起 始 密码 子 、 长 度 大 于 300 pb 的 非 
重复 序列 进行 中 性 绘图 分 析 结 果 显 示 (图 3) ， 宁 冈 青冈 、 曼 青冈 、 毛 曼 青 冈 、 竹 叶 青冈 
GC3 和 GCio 的 相关 系数 分 别 为 0.007 93、0.053 84、0.005 495、0.06， 回 归 系 数 分 别 为 0.116 
33. 0.336 92、0.342 67. 0.356 94， 双 尾 检验 结果 均 为 P<0.01， 显 示 出 GCio 和 GC3 之 间 的 
相关 性 较 弱 ， 组 成 密码 子 碱 基 的 第 1、2 位 和 第 3 位 之 间 存 在 显著 差异 。 推 测 自然 选择 对 4 
种 青冈 亚 属 植物 叶绿体 基因 组 密码 子 的 使 用 偏好 性 的 影响 更 大 。 


密码 子 第 1、2 位 平均 GC 含量 
on the ls 


Averag 


70.116 33x+0.395 32 
R=0.007 93 


pa 
M 0.1 0.2 0.3 04 0.5 0.6 


密码 子 第 3 位 GC 含量 
GC content on the 3rd site (GC3) 


70.342 67x+0.333 89 
R=0.054 95 


M 0.1 02 0.3 04 0.5 0.6 


密码 子 第 3 位 GC 含量 
GC content on the 3rd site (GC3) 


内 叶绿体 基 


A. 宁 冈 青冈 叶绿体 基因 组 中 性 绘图 ; B. 曼 
图 ; D. 竹叶青 冈 叶 绿 体 基因 组 中 性 绘图 。 


y70.336 92x+0.335 65 
R7-0.05384 


P d 


MY 0.1 02 0.3 04 0.5 0.6 


密码 子 第 3 位 GC 含量 
GC content on the 3rd site (GC3) 


3=0.356 94x+0.329 79 
R=0.06 


MT] 0.1 02 03 04 0.5 0.6 


密码 子 第 3 位 GC 含 量 
GC content on the 3rd site (GC3) 


[KH H 


PERR; C. 毛 曼 青冈 叶绿体 基因 组 中 性 绘 


A. Neutrality polt analysis of the chloroplast genome of Quercus ningangensis; B. Neutrality polt analysis of the 


chloroplast genome of Q. oxyodon; C. Neutrality polt analysis of the chloroplast genome of Q. gambleana; D. 


Neutrality polt analysis of the chloroplast genome of Q. neglecta. 
图 3 中 性 绘图 分 析 
Fig. 3 Analysis of neutrality plot 


2.2.3 ENC-plot £2 


4 种 青 闷 亚 属 植 物 的 ENC-plot 分 析 结果 显示 《图 


准 曲 线 之 下 , 说 明 


分 析 


自然 选择 对 这 部 分 基因 的 密码 子 使 用 偏好 性 的 影响 更 大 ; 也 有 少 部 分 基因 


4) ， 大 部 分 基因 的 实际 ENC 值 在 标 


的 实际 ENC 值 在 标准 曲线 附近 或 上 方 ， 说 明 突 变 对 其 密码 子 的 使 用 偏好 性 影响 更 大 。ENC 
比值 频率 分 布 ( 表 5) 显示 ， 宁 内 青冈 、 曼 青冈 、 毛 曼 青 网 、 竹 叶 青 内 中 位 于 -0.05~0.05 区 


间 的 基因 分 别 有 19、18、18、19 个 ， 分 别 占 37%、35%、35%、37%， 
因 的 影响 较 大 ,在 这 4 种 青冈 】 


表明 突变 对 这 些 基 


严 属 植物 中 还 有 33-34 个 基因 的 实际 ENC 值 与 理论 ENC 值 差 


距 较 大 ， 这 部 分 基因 受 自然 选择 的 影响 作用 较 大 。ENC-plot 分 析 结 果 和 ENC 比值 频率 分 布 
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显示 , 4 种 青冈 亚 属 植物 的 叶绿体 基因 组 密码 子 偏好 性 同时 受 突变 和 自然 选择 的 影响 ,但 总 
的 来 说 ， 自 然 选 择 的 影响 大 于 突变 带 来 的 影响 。 


g g 
á A 
E z o5 
s$ ui 
Es gs 40 
LE. XH 
EE EE 
Ea ma 3 
B o2 FEET 
Ei Es 
10 10 
Oo or 02 03 04 05 06 07 08 09 T0 Oo or 02 03 04 os 06 07 08 09 TO 
密码 子 第 3 位 GC 含量 密码 子 第 3 位 GC 含量 
GC content on the 3rd site (GC3) GC content on the 3rd site (GC3) 
70. 70. 
e D 
- 2060 
加 g 
加 á 
E 2 so 
FE EE 
T3 ES 
ES z $ 40 
€i $5 
E F 
8$ ES 30 
goo 8 om 
Es Es 
10 10 
Wo wi 02 03 04 03 06 07 08 09 To Oo or 02 03 04 03 06 07 08 09 10 
密码 子 第 3 位 GC 含量 密码 子 第 3 位 GC 含量 
GC content on the 3rd site (GC3) GC content on the 3rd site (GC3) 


A. 宁 冈 青冈 叶绿体 基因 组 ENC-plot 分 析 ; B. 曼 青 冈 叶 绿 体 基因 组 ENC-plot 分 析 ; C. 毛 曼 青冈 叶绿体 基 
因 组 ENC-plot 分 析 ; D. 竹叶青 冈 叶绿体 基因 组 ENC-plot 分 析 。 
A. ENC plot analysis of the chloroplast genome of Quercus ningangensis; B. ENC plot analysis of the chloroplast 


PI 


genome of Q. oxyodon; C. ENC plot analysis of the chloroplast genome of Q. gambleana; D. ENC plot analysis of 
the chloroplast genome of Q. neglecta. 
4 ENC-plot 分 析 
Fig. 4 Analysis of ENC-plot 


表 5ENC 比值 频数 分 布 
Fig. 4 Distribution of ENC ration 


物种 组 限 组 中 值 频数 组 频 
Species Class boundary Group mid-value Frequency number Frequency 
-0.15—-0.05 -0.08 3 0.06 
TP -0.05—0.05 0.01 19 0.37 
Quercus ningangensis 0.05—0.15 0.06 19 0.37 


0.15—0.25 0.19 10 0.19 


zT 


Q. oxyodon 


毛 曼 青冈 
Q. gambleana 


竹叶青 冈 
Q. neglecta 


2.2.4 PR2-plot 绘图 分 析 


PR2-plot 绘图 中 (图 5 


中 心 的 点 则 受 自然 选择 的 影响 。4 种 青 
均匀 ,大 部 分 基因 分 布 于 右 下 4 


0.25—0.35 
-0.15—-0.05 
-0.05—0.05 
0.05—0.15 
0.15-0.25 
0.25—0.35 
-0.15—-0.05 
-0.05—0.05 
0.05—0.15 
0.15-0.25 
0.25—0.35 
-0.15—-0.05 
-0.05—0.05 
0.05—0.15 
0.15—0.25 
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突变 的 影响 ， 大 多 远离 


内 亚 属 植物 的 PR2-plot 绘图 分 


析 显 示 ， 基 因 的 分 布 不 
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J, WHH 4 种 青冈 亚 属 植物 叶绿体 基因 组 中 密码 子 第 3 位 碱 基 
表明 4 种 青冈 亚 属 植物 密码 子 的 偏好 性 受到 了 突变 和 自然 选择 
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A3. T3. G3. C3 分 别 表 示 密 码 子 第 三 位 上 各 个 碱 基 的 含量 。A. 宁 冈 青冈 叶绿体 基因 组 PR2-plot 分 
Wr; B. 曼 青 冈 叶绿体 基因 组 PR2-plot 分 析 ; C. 毛 曼 青冈 叶绿体 基因 组 PR2-plot 分 析 ; D. 
竹叶青 内 叶绿体 基因 组 PR2-plot 分 析 。 
A3, T3, G3, and C3 indicate the content of each base in the 3rd position of the codon. A. PR2-plot of chloroplast 


genes in Quercus ningangensis; B. PR2-plot of chloroplast genes in Q. oxyodon; C. PR2-plot of chloroplast genes 


in Q. gambleana; D. PR2-plot of chloroplast genes in Q. neglecta. 


图 


5 PR2-plot 分 析 


Fig. 5 Analysis of PR2-plot 


202208.00040v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


2.3 叶绿体 基因 组 差异 比较 分 析 

选 宁 冈 青冈 的 叶绿体 基因 组 做 参考 ， 利 用 mVISTA 对 曼 青 网 、 毛 曼 青冈 、 竹 叶 青 冈 、 
小 叶 青 冈 、 福 叶 青冈 、 西 畴 青冈、 青冈 的 叶绿体 基因 组 序列 进行 可 视 化 比较 分 析 。 从 分 析 结 
果 可 以 看 出 青冈 亚 属 植物 的 叶绿体 基因 组 具有 较 高 的 保守 性 , 主要 的 变异 区 域 位 于 非 编 码 区 ， 
而 ycfl 基因 在 蛋白 质 编码 区 域 发 生变 异 较 大 (图 6) 。LSC 区 域 的 变异 较 大 ，IRs 区 域 变异 
小 较为 保守 ，SSC 区 域 的 变异 小 于 LSC 区 域 。8 种 青冈 亚 属 植物 的 rRNA 基因 高 度 保守 ， 
几乎 没有 发 生变 异 (图 6) 。 
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64 种 青冈 亚 属 植物 及 近 缘 种 叶绿体 基因 组 的 可 视 化 比 对 


Fig. 6 Visualization alignment of four chloroplast genomes in 


Quercus subg. Cyclobalanopsis and its related species 

24 IR 边界 分 析 

使 用 在 线 软件 IRscope 分 析 4 种 青冈 亚 属 植物 及 其 近 缘 种 叶绿体 基因 组 的 TR. 扩张 和 收 
缩 情况 ， 以 及 识别 LSC、IR、SSC 的 边界 位 置 (图 7) 。 在 8 种 青冈 中 ， 各 边界 附近 分 布 的 
基因 一 致 ， 但 距 边界 的 距离 略 有 不 同 。 其 中 ，7ps79 基因 距 SSC/IRb 边界 均 为 11 bp， 西 畴 青 
内 和 青冈 的 tend. 基因 距离 IRa/LSC 边界 的 距离 为 16 bp, 其 余 青 冈 为 1bp。 在 毛 曼 青冈 和 竹 
叶 青 四 中 ndhF 基因 ， 有 一 个 碱 基 位 于 Rb 区 域 ， 其 余 2255 bp 均 位 于 SSC KR, HRE 
中 ndhF 基因 的 长 度 为 2255 bp， 且 均 位 于 SSC 区 域 。 高 等 植物 的 叶绿体 基因 组 中 ycf1 基因 
为 第 二 长 的 基因 ， 其 横 跨 这 8 种 青 闷 叶 绿 体 基因 组 的 SSC/IRa 边界 ,在 宁 冈 青冈 、 小 叶 青冈 、 
褐 叶 青冈 、 西 畴 青冈 、 青 冈 中 ，)xcej 基因 有 1 060 bp 位 于 Ra 区 域 ， 在 曼 青 冈 、 毛 曼 青 冈 、 
竹叶青 冈 中 分 别 有 1 056 bp. 1054 bp. 1062 bp 位 于 Ra Ki, 在 8 种 青冈 的 IRb 区 域 也 分 
别 拷贝 了 相应 长 度 的 yej 片段 ， 这 部 分 被 视 为 ycf7 假 基 因 Coycfl) . 8 种 青冈 中 ， vycfl 
基因 横 跨 IRb/SSC， 有 56-58 bp 位 于 IRb/SSC 上 游 。 
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Fig. 7 Comparison of LSC, SSC, and IR border regions among four chloroplast genomes in 


Quercus subg. Cyclobalanopsis and its sibling species 


2.5 系统 发 育 分 析 


选取 NCBI 上 公开 释放 的 8 种 壳 斗 科 植 物 的 叶绿体 基因 组 与 该 研究 组 装 的 4 种 青冈 亚 属 
植物 的 叶绿体 基因 组 序列 ， 以 台湾 水 青冈 Fagus hayatae) 和 米 心 水 青 风 CF engleriana) 
作为 外 类 群 ， 分 别 用 ML 法 和 NJ 法 构建 系统 发 育 树 ， 来 探讨 青冈 亚 属 和 栎 亚 属 之 间 的 关系 
及 毛 曼 青 由 和 曼 青 冈 的 关系 ,并 为 青冈 亚 属 系统 发 育 问题 的 解决 提供 基础 资料 。 两 种 方法 构 


建 的 系统 发 育 树 均 得 到 了 较 高 的 支持 率 且 获得 的 拓扑 结构 一 致 。 系 统 发 育 树 显示 (图 8) ， 


水 青冈 属 (Fagus) 和 三 棱 栎 属 CFormanodendron) 位 于 系统 发 育 树 的 基部 ， 为 克 斗 科 中 较 
早 分 化 出 来 的 类 群 。 曼 青冈 和 毛 曼 青冈 成 姐妹 群 关 系 ， 毛 曼 青 内 为 4 种 青冈 亚 属 植物 中 最 早 
分 化 出 来 的 类 群 。 栎 亚 属 并 未 形成 一 个 单 系 群 ， 川 溃 高 山 标 〈Q. aquifolioides) 与 青冈 亚 属 


ChinaXiv 合 作 期 刊 


87/99 Q SEP Quercus oxyodon ON258628 
8897| L @ 竹叶青 网 O ON258631 
Q FARER Q. ON303301 青冈 亚 属 
Q EZEN Q. ON258629 Quercus subg. Cyclobalanopsis 
倒 卵 叶 青冈 O. obovatifolia MG356785.1 
Ju E RS IUBE O. aquifolioides KX911971.1 
刺 叶 高 山 栎 O. spinosa MG678038.1 栎 亚 属 
EHR Q. engleriana MZ196209.1 
FEER Q. variabilis NC 031356.1 


99/100 
100/100 


100/100 


100/100 


subg. Quercus 


ZIRH Trigonobalanus doichangensis NC 023959.1 
BAKE] Fagus hayatae MW846258.1 
米 心 水 青冈 F. engleriana NC 036929.1 


ik: 上 面 的 节点 数 是 支持 值 ， 左 边 是 ML 自 展 值 ， 是 NJ 自 展 值 。 


Note: Number above nodes are Support values with ML bootstrap values on the left, NJ bootstrap values on the 


right. 


图 8 基于 ML 法 和 NJ 法 对 12 个 物种 叶绿体 基因 组 序列 构建 系统 发 育 树 
Fig. 8 Phylogenetic relationships of the 12 species inferred from ML and NJ methods 


3.4 叶绿体 基因 组 比较 及 密码 子 偏好 性 分 析 

被 子 植 物 的 叶绿体 全 基因 组 长 度 在 120 000~170 000 bp 之 间 (Tangphatsornruang et al., 
2010), 本 研究 中 宁 冈 青冈 、 曼 青冈 、 毛 曼 青 冈 、 竹 叶 青 冈 叶绿体 基因 组 的 长 度 为 160 784~160 
906 bp， 在 叶绿体 基因 组 的 大 小 、 结 构 、 基 因 的 数量 、 总 的 GC 含量 及 4 个 区 域 的 长 度 和 
GC 含量 上 都 高 度 保守 。 青 冈 亚 属 的 叶绿体 基因 组 基本 特征 与 其 他 壳 斗 科 植 物 相 似 ， 壳 斗 科 
与 其 他 科 物 种 相 比 基 因 密度 更 大 ， 对 环境 的 适应 更 强 ， 稳 定性 更 高 〈 高 潇潇 , 2020) 。 从 4 
种 青冈 亚 属 植物 的 叶绿体 全 基因 组 比较 分 析 可 以 看 出 ,4 种 青冈 亚 属 植物 的 叶绿体 基因 组 高 
度 保守 ， 变异 区 域 主要 在 非 编码 区 。 青 冈 亚 属 植物 的 鉴定 主要 依据 壳 斗 的 特征 ， 在 没有 果实 
的 情况 下 ， 部 分 种 难以 鉴定 ， 可 以 利用 叶绿体 基因 组 筛选 变异 较 大 的 区 域 ， 开 发 成 DNA 条 
形 码 ， 便 于 青 闷 亚 属 植物 的 鉴定 。 通 过 对 8 种 青冈 亚 属 植 物 叶 绿 体 基因 组 的 研究 ， 篇 选 出 
trnH-GUG~trnQ-UUG , ropClI-petN. trnT-UGC~trnM-CAU, ycf4~psbJ 及 ycfl 区 域 作为 DNA 
条 形 码 ， 为 后 续 青 闷 亚 属 植物 的 鉴定 提供 便利 。 被 子 植物 叶绿体 基因 组 IR. 区 域 的 扩张 和 收 
缩 普遍 存在 于 植物 进化 过 程 中 (Hansen et al., 2007; Davis & Soreng, 2010; Huang et al., 2014)。 
青冈 亚 属 植物 位 于 IR 区 域 的 基因 相同 ,说 明 青 闷 亚 属 IR. 边界 较为 保守 。4 种 青冈 亚 属 植物 
JEI yef1 假 基 因 (了 ycf1〉 ， 在 壳 斗 科 标 亚 属 ， 栗 属 〈Castanrea) 部 分 种 中 也 出 现 了 ycfl 
基因 在 IRb/SSC 拷贝 为 假 基 因 的 现象 ( 杨 颜 慈 , 2018; 高 潇潇 , 2020) . 

密码 子 使 用 偏好 性 是 生物 进化 的 重要 现象 , 主要 受 基因 组 的 大 小 、 碱 基 突 变 、 基因 漂 变 、 
自然 选择 、 基 因 的 表达 水 平 、tRNA 丰 度 、GC 含量 、 蛋 白质 结构 等 因素 影响 (Romero etal., 
2000; Duret , 2000; Angellotti et al., 2007) 。 从 GC3 的 含量 可 以 看 出 本 研究 中 的 青冈 亚 属 
植物 更 偏好 以 AT 结尾 的 密码 子 ， 这 一 现象 在 被 子 植 物 中 广泛 存在 (Clegg et al., 1994; 


x 


Tangphatsornruang et al., 2010; Delannoy et al., 2011) 。 同 义 密 码 子 相对 使 用 度 结果 显示 ， 青 
闷 亚 属 和 栎 亚 属 植物 的 同 义 密码 子 使 用 度 相 似 ， 米 档 、 长 果 柯 中 UAA 的 使 用 频率 较 所 列 旧 
其 他 植物 低 , 说 明 壳 斗 科 植物 的 密码 子 使 用 偏好 性 相似 , 但 不 同 的 属 之 间 仍 存在 着 一 定 的 差 
异 。3 种 绘图 分 析 结 果 显 示 , 4 种 青 闷 亚 属 植物 的 密码 子 偏好 性 受 选择 压力 和 突变 同时 影响 ， 
但 选择 压力 为 主要 影响 因素 ， 这 与 米 档 CCastanopsis carlesii) 〈 江 淑 珍 等 , 2021) 、 降 香 黄 
T8 (Dalbergia odorifera) (JAIR IESE, 2021), IT RHE (Erigeron breviscapus) (Æ TKE, 2021) 
等 得 出 的 结果 相似 。 
32 系统 发 育 关 系 及 种 间 关 系 分 析 
系统 发 育 树 结果 显示 ， 水 青冈 属 (Fagus) MEZIRE Clrigonobalanus) 分 化 较 早 ， 与 

前 人 根据 化 石和 地 理 信息 系统 对 充 斗 科 起 源 等 研究 的 结果 一 致 〈 周 浙 昆 , 1999) 。 前 人 根据 
花粉 特征 及 核 基因 将 川 汗 高 山 标 、 刺 叶 高 山 标 、 巴 东 栎 归 为 lex 组 、 栓 皮 栎 归 为 Cerris 组 
(Hubert F et al., 2014) ， 但 核 标 记 、 质 体 标 记 的 研究 结果 中 显示 ，Ilex 组 的 亚洲 种 常 出 现 
穿插 在 Cerris 组 中 的 现象 (Simeone MC et al., 2013; Hubert F et al., 2014) ， 本 研究 中 ，4 
种 栎 亚 属 植物 并 未 形成 一 个 单 系 群 , 刺 叶 高 山 栎 (Quercus spinosa) 5 ERER CO. engleriana), 
FEER CO. variabilis) 聚 在 一 枝 ，Hlex 组 和 Cerris 组 之 间 出 现 穿插 的 现象 ， 与 核 标 记 、 质 体 
标记 的 结果 一 致 。 在 《中 国 植物 志 》《 陈 焕 镭 和 黄 成 就 , 19980. 、《 中 国 树木 志 》《〈 郑 万 钩 ， 
1985) 、《 云 南 植物 志 》 〈 中 国 科 学 院 昆明 植物 研究 所 , 1979) 、《 贵 州 植物 志 》 CEN 
和 李 永 康 , 1982) 、 等 地 方 植物 志 中 ， 毛 曼 青 由 和 曼 青 冈 是 作为 两 个 单独 的 种 处 理 。 罗 艳 和 
周 浙 昆 (2001a) 、 邓 敏 〈2007) 等 根据 曼 青冈、 毛 曼 青 冈 叶 背 的 星 状 毛 等 特征 ， 认 为 毛 曼 
青冈 应 作为 曼 青 冈 的 一 个 变种 处 理 。 而 刘 凌 燕 等 2008) 根据 聚 类 分 析 的 结果 认为 应 将 毛 曼 
青冈 作为 长 叶 青 冈 〈O. longifolia) 的 变种 处 理 。 本 研究 中 系统 发 育 分 析 结 果 显 示 ， 曼 青 囚 
首 叶 青冈 亲缘 关系 较 近 , 曼 青 风 和 毛 曼 青冈 未 形成 姐妹 群 。 毛 曼 青 由 和 曼 青 内 的 果实 和 叶 
HAH, 叶 背 毛 的 颜色 和 形态 存在 差异 , 但 在 不 同 地 区 采集 的 叶 背 毛 差距 很 大 ( 邓 敏 , 2007)。 
扫描 电镜 下 观察 两 者 的 花粉 外 壁 纹 饰 纹 饰 上 也 存在 区 别 , 毛 曼 青冈 的 花粉 外 壁 纹饰 为 颗粒 
天 或 聚合 疣 状 ， 曼 青冈 的 花粉 外 壁 为 颗粒 状 纹饰 《 王 萍 莉 和 溥 发 易 , 2004) 。 本 研究 中 采 
的 毛 曼 青冈 叶 密 被 黄 褐色 星 状 毛 ， 而 曼 青 冈 叶 背 的 星 状 毛 为 白色 且 较 毛 曼 青 内 稀 疏 。 综 上 
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KIRE FESE PURI E IBAMA. EAE HIE RUREL UE 
曼 青 风采 自 广西 ， 两 地 相距 较 远 , 我 们 推测 可 能 人 
流 受 阻 等 因素 导致 本 研究 中 的 毛 曼 青 网 和 曼 青 冈 逐 渐 形 成 两 个 种 。 
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